Percolation in random environment 
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We consider bond percolation on the square lattice with perfectly correlated random probabilities. 
According to scaling considerations, mapping to a random walk problem and the results of Monte 
Carlo simulations the critical behavior of the system with varying degree of disorder is governed 
by new, random fixed points with anisotropic scaling properties. For weaker disorder both the 
magnetization and the anisotropy exponents are non-universal, whereas for strong enough disorder 
the system scales into an infinite randomness fixed point in which the critical exponents are exactly 
known. 



I. INTRODUCTION 

Percolation is a paradigm for random processes , in 
which the i-th bond (or site) of a regular lattice is oc- 
cupied with a probability, pi, which is generally taken 
independent of its position, pi — p. In percolation theory 
one is interested in the properties of clusters, in partic- 
ular in the vicinity of the percolation transition point 
p = pc, when clusters with diverging size are formed. 
Using a close analogy with thermal phase transitions, 
which is based on the Q — > 1 limit of the ferromagnetic 
Q-state Potts model [|| , a scaling theory has been devel- 
oped and in two dimensions many, conjecturedly exact 
results have been obtained by conformal field theory ||] 
and by Coulomb-gas methods [Q. 

In real systems, however, the occupation probabilities 
are generally inhomogeneous, i.e., position, direction or 
neighborhood dependent, and there are some correlations 
between them. The effect of quenched disorder, i.e. when 
the occupation probabilities are position dependent ran- 
dom variables, can be studied by scaling considerations. 
According to the Harris criterion Q the relevance or ir- 
relevance of the effect of quenched disorder on the perco- 
lation transition depends on the sign of the specific heat 
exponent, a, of the corresponding pure Potts model in 
the Q — > 1 limit. Since in any dimension a < Q, the 
critical properties of ordinary and "random" percolation 
are equivalent. Another form of perturbations, e.g. long- 
range correlations between occupation probabilities |^ or 
anisotropy, such as in directed percolation however, 
leads to modified critical properties. 

In the present paper we consider the combined effect 
of disorder, anisotropy and correlations, when the occu- 
pation probabilities are random variables, which are per- 
fectly correlated in a dd dimensional subspace. This type 
of behavior could be relevant to describe the properties 
of oil or gas inside porous rocks in oil reservoirs, when 
the rock has a layered structure 

Models with perfectly correlated disorder play an im- 
portant role in statistical physics and in the theory of 
(quantum) phase transitions. Among the early work we 



mention the partially exact solution of the McCoy- Wu 
model (which is the two-dimensional Ising model with 
layered randomness) and the field-theoretical investiga- 
tions by Boyanovski and Cardy |^ . As a matter of fact in 
random quantum systems disorder is perfectly correlated 
along the (imaginary) time direction, i.e. here = 1. 
For these systems, in particular for random quantum 
spin chains, i.e. in (1 + 1) dimension, many new, pre- 
sumably exact results have been obtained recently by a 
strong disorder renormalization group (SDRG) method 
JlOt . It was found that for strong enough initial disorder 
the critical behavior of several systems is governed by 
a so-called infinite randomness fixed point (IRFP) [ pd] |, 
with unusual scaling properties. Here we mention recent 
calculations on the random transverse-field Ising model 
(RTIM) |T2|jr^ ], random quantum Potts and clock mod- 
els [14], random antiferromagnetic Heisenberg spin chains 
p^,16| and ladders and also non-equilibrium phase 
transitions in the presence of quenched disorder |l8| . In 
many cases a cross-over between weak and strong dis- 
order regimes has been observed and a general scaling 
scenario has been proposed . 

In the present paper we study percolation on a square 
lattice with strip random occupation probabilities. We 
investigate the critical behavior of the system with vary- 
ing strength of disorder by scaling considerations, by ran- 
dom walk mappings and by Monte Carlo (MC) simula- 
tions. The structure of the paper is the following. The 
model and the relevant physical quantities are introduced 
in Sec. II. Investigations in the weak and strong disorder 
limits are given in Sec. HI., MC simulations for inter- 
mediate disorder are presented in Sec. IV. The paper is 
closed by a discussion in Sec. V. 



II. THE MODEL 

We consider bond percolation on a square lattice with 
sites {j, j}, 1 < * < i and 1 < j < K , where the oc- 
cupation probabilities, < p < 1, are random variables, 
which are perfectly correlated along vertical lines, as 
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indicated in Fig. |^. If the average value of the occupa- 
tion probabihties exceeds a critical value, (p) > Pc^ there 
is a percolation transition in the system. The value of 
Pc can be determined by noticing that under a duality 
transformation, which maps the ordered and the disor- 
dered phases of the system into each other, the layered 
structure of the system is preserved and the dual value 
of the local probability is transformed as: p = 1 — p 
pOf . Consequently the probabihty distribution, P{p), is 
transformed into P{p) = P(l— p) and the random system 
is self-dual, if the probability distribution is symmetric: 
P{p) = P{1 — p) and thus the average value of p is given 
by (p) = Pc = 1/2. Since there is one phase transition 
in the system, the self-duality point corresponds to the 
critical point and the distance of the critical point, t is 
defined as: 



i = (p) " Pc 



(1) 



K 



2 



1 = 



1 



L 



FIG. 1. Percolation on the square lattice with random 
bond occupation probabilities, which are perfectly correlated 
along the vertical direction. The portion of the lattice hav- 
ing the same occupation probability, p, is denoted by bold 
hnes, whereas the corresponding part of the dual lattice with 
p = 1 — p is shown by dashed lines. 

In the presence of quenched disorder the mean value of 
a physical observable, <&, is calculated as [($)]av, where 
(. . .) denotes thermal averaging for a given realization of 
the disorder and [. . .]av stands for disorder averaging. In 
percolation the basic quantities of interest are the frac- 
tal and connectivity properties of the largest clusters. In 
the following we use the concept of anisotropic scaling 
when the correlation lengths in the two directions, 
(which correspond to the extensions of the largest clus- 
ters) involve different critical exponents: ~ 1^1 ""^"^ ^-nd 
(Cii ^ i"'^" ■ Thus the anisotropy exponent 



(2) 



is generally different from one. In the ordered phase, 
i > 0, the number of points belonging to the infinite 
cluster, A'o, scales around the transition point as: 



(3) 



where [3 is the critical exponent of the order parameter. 
At the critical point, t = 0, fixing the ratio KjU' = 0{1) 



we obtain: 



(4) 



where the two fractal dimensions of the infinite cluster 
are given by: 



d_L = 1 + z — P/i'±_ 



(5) 



and = d^/z. The distribution of cluster sizes, R(N), 
at the critical point asymptotically behaves as: 



R{N)dN = N-''R{N/L'^^)dN , 



(6) 



where r = 2 + (3/{i^±d±). This relation can be obtained 
by generalizing the similar result for ordinary percolation 

0. 

Correlation between two sites with coordinates, {ii,ji} 
and {z2, j2}, is defined as the expectation value of the 
connectivity, (^Cfii, ji}, {j2, J2}), which is 1, if the two 
sites belong to the same cluster and zero otherwise. Here 
we mainly consider correlations in the perpendicular di- 
rection 



|av 



(7) 



where an average over the vertical coordinate, ji — j2 — j 
is also performed. When correlations in the bulk are cal- 
culated we use periodic boundary conditions (b.c), (thus 
i — L + 1 = 1), take maximal distance between the sites, 
^2 = ii + L/2, and average over the position ii. The 
average bulk correlations, calculated in this way, scale at 
the critical point as: 



CiiL) 



(8) 



where r]± = 2(3 /v^_. We also considered the system with 
free boundaries at i = 1 and i = L and calculated the 
correlations between two surface sites. This end-to-end 
correlation function at the critical point asymptotically 
behaves as: 



C^il,L) = Cl{L)^ L- 



(9) 



where the decay exponent, rj'^, is related to the surface 
fractal properties of the infinite cluster. Closing this sec- 
tion we quote the values of the critical exponents for two- 
dimensional ordinary percolation [Q: 



4/3, 77^°) = 5/24, 77 
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(10) 
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III. STRENGTH OF DISORDER: LIMITING 

CASES 

The strength of disorder, A, is related to the broad- 
ness of the probabiUty distribution, P{p). In terms of the 
integrated probabihty distribution, n(p) = P{p')dp' , 
we introduce the probabihties, Pi/4 and P3/4 with the 
definitions: n(pi/4) = 1/4 and 11(^3/4) = 3/4. Since the 
central half of the distribution is located in the region: 
Pi/4 P l£ P3/4 its relative width is measured by: 



A = 



P3/4 ~ Pl/4 
1 - P3/4 + -Pl/4 



(11) 



what we can identify with the strength of disorder. 

In this paper we used two specific forms of the distri- 
bution. For the bimodal distribution (0 < q < 1/2, q — 
1-9): 



Pbin [p] 



1 -t 



5{p-q) 



1 + t 



5{p-q) 



(12) 



the critical point is located aX t — Q and the strength 
of disorder is given by Abi„ = (1 — 2q)/2q. Thus, as 
expected the bimodal disorder is weak for 5 « 1/2 and 
strong for g ^ 1/2. 

The other distribution we use has a power-law form: 



Ppow {p^ 



1 



D2p \p 



-1+1/D 



0<p<p<l, (13) 



and Ppow{l 



pow{p)p/{l - p), for p < p < 1. 



P) ^ Pj- 

The distance from the critical point is measured by 
t = {p- l/2)/{D + 1), and for p = 1/2, i.e. for i = 
the distribution is indeed symmetric. In this case the 
strength of disorder is given by: Apow = 2^ — 1. Thus 
for D = we recover the ordinary percolation and the 
strength of disorder is monotonically increasing with D. 
Therefore D will be often called as the disorder parame- 
ter of the distribution. 



A. Weak disorder 

In the limit of weak disorder one usually decides about 
the relevance- irrelevance of the perturbation by perform- 
ing a stability analysis at the ordinary percolation fixed 
point. Generalizing the method by Harris j|] the cross- 
over exponent due to correlated disorder is calculated as: 



.(0) 



2/3 



(14) 



where we used i^^^-' = 4/3 in Eq.([T^). Since > 0, even 
weak correlated disorder is a relevant perturbation, thus 
a new random fixed point is expected to control the crit- 
ical behavior of the model. 



B. Strong disorder: Mapping to random walks 

Next we turn to study the behavior of the system for 
extremely strong disorder using the bimodal distribution 
in Eq.(p^ in the limit g — s- 0. In this limiting case the 
percolation in a given layer with a probability pi has a 
simple, anisotropic structure (for an illustration see Fig. 
H). If this probability is extremely large, pi = q, then 
here almost all bonds are occupied, except of a very small 
fraction of q. Since the typical distance between two non- 
occupied bonds is I ~ 1/q, the cluster in the z-th layer is 
composed of long connected units of typical size I. On the 
other hand, if the probability is extremely small, pj — q, 
then almost all bonds in this layer are unoccupied, except 
of a very small fraction of q. Since the typical distance 
between two occupied bonds is I 1/q, the cluster in 
the j-th column is composed from long empty units of 
typical size I. Notice the duality in the structure of the 
two types of column. 



n, 



q q q q q 

FIG. 2. Structure of the percolation cluster in the extreme 
bimodal distribution with pi = P2 = q and ps = P4 = Ps = 5 
(here with q « 1/4). In a layer with extremely large 
(small) probability there are connected (empty) units of typ- 
ical length I ~ 1/q. The number of sites of the connected 
cluster at the other surface of a strip of width, k, Uk is given 
by: ni ~ 1/q, n2 ~ 1/q'^, ns ^ 1/q and n4 = 0(1) (see text). 
In the limit q — > the cluster ends at = 4, thus ns = 0. 

With this prerequisite we consider the order-parameter 
in the surface column, (L) , which is the fraction of sur- 
face sites belonging to a cluster of horizontal extent L. 
In order to make a statement about the value of ms{L) 
we consider parallel strips of width k < L and introduce 
the quantity, Uk, as the typical number of bonds at the 
fc-th (i.e. surface) column of a cluster, which is connected 
to the other surface of the strip. Starting with fc = 1 we 
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have two possibiHties. For extremely small probability, 
Pi = q, there is no surface cluster in the system, thus we 
have III = 0. Otherwise, for pi = q, a surface site is con- 
nected to all sites of a "connected unit" of length I, thus 
we have ni I 1/ q. For fc = 2, if the probability is 
extremely large in the second layer, too, pi = P2 = q, 
then a surface cluster extends up to the second layer 
and its vertical size, which is given by n2, can be esti- 
mated as follows (see Fig. |2|). The end of a cluster is 
signalled by the fact that in both columns unoccupied 
bonds are in neighboring positions, which happens with 
a probability q^ , from which the typical size of a cluster 
n2 ~ l/f?^ ~ TT-i/q, follows. Repeating this argument for 
Pi = g, i = 1, 2, . . . , /c we obtain Uk ^ l/g*'' ~ Uk-i/q. 
Now having a small probability at the following layer. 
Pi = q, i = 1,2, ... ,k and Pk+i — q, then only a fraction 
of q of the sites have a further connection, thus will 
be reduced by a factor q giving rik+i ^ n^q. Inclusion 
of any further layer with an extremely small probabil- 
ity will reduce rij by a factor of q, until we arrive at 
rij' — 0{1), when for the next small probability layer 
we have n^'+i = 0, thus the surface cluster ends at this 
distance. 

From this example we can read that the Uk num- 
bers are either integer powers of rik ~ 1/q-^'', for 
Xk = 0, 1, . . ., or Tifc = 0, if formally Xk < 0. Further- 
more, we have the transformation rules: 



nk+i 



nk/q, Pk+i = q 
rikq, Pk+1 = q 



(15) 



where in the second case Uk+i = 0, if Uk = 0(1). At this 
point we can formulate the condition that the surface 
magnetization in a given sample (in a rare realization) is 
m^(i) 0(1), if Uk > 0(1), for all fc = 1, 2, . . . L. For 
all other cases ms{L) = 0. Consequently to calculate the 
average value of ms{L) it is enough to find the fraction 
of rare realizations, for which ms{L) — 0(1), since 
[jTisJav Pl- To calculate p'l we use a random walk 
(RW) mapping (see an illustration in Fig. in which 
to each disorder realization we assign a one-dimensional 
RW, which starts at Xq — and takes its fc-th step 
upwards, Xk — 1 (downwards, Xfc = — 1 ) if the corre- 
sponding bond occupation probability is extremely large, 
Pk — q (extremely small, pk — q). The position of the 
walker at the fc-th step, Xk — J2i=i is related to Uk as 
Uk ~ q~^'' . Then, as argued before, the surface cluster 
extends up to a vertical distance, L, if Xk > 0, for every 
k = 1, 2, ... L, i.e. the RW has a surviving character. 

At the critical point of the percolation problem, t — 0, 
the corresponding RW is unbiased, and the fraction of 
surviving L-step RW-s scales as ~ L~^/^. Now the 
fraction of clusters which connect the two free boundaries 
of the strip over a distance L, and thus contribute to 
the average end-to-end correlations in Eq. (^ , is given by 
since at each site there should be an independent 



percolating surface cluster, which meet in the middle of 
the system. Consequently the average end-to-end corre- 
lations at the critical point scale as C^{L) ~ L^^ thus 
the corresponding decay exponent in the strong disorder 
limit is given by: 



s.(oo) 



1 



(16) 




k= 1 



FIG. 3. Illustration of the RW mapping of percolation for 
a given realization of the extreme binary distribution. Layers 
with high, q, (low, q) probability are drawn by thick (thin) 
lines and the corresponding RW makes a step of unit length 
upwards (downwards). The position of the RW, in the fc-th 
step, Xk, is related to, Uk, the number of typical sites in the 
fc-th layer of percolation, which are connected to a given sur- 
face site as Uk ~ The surface cluster extends to a 
distance, L, if Xk > 0, for all fc = 1, 2, . . . L, thus the RW has 
a surviving character. 

Using the RW mapping one can easily estimate the per- 
pendicular size of the percolating clusters, which is given 
by ^ ^ q^'^/'^. Since the transverse fluctua- 

tions of unbiased surviving RW-s scale as ^ L^^"^ 

we obtain in the strong disorder limit 



(17) 



in Eq. 



IS 



Consequently the anisotropy exponent, z 
formally infinite for strong disorder. 

Another results can be simply obtained by noticing 
that the same type of RW mapping applies to the one- 
dimensional RTIM too, so that we can simply bor- 
row the results obtained in this case. 

For bulk correlations one should consider the fraction 
of realizations, pL, for which a given bulk site belongs 
to a connected cluster of vertical size, L. As was shown 
in for these realizations the thermal average of the 
position of the RW has a surviving character. The frac- 
tion of these walks is given by |p^,p3l: pl ~ L-(3-v^)/4^ 
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consequently the critical average bulk correlations being 
C^{L) ~ (pl)^ have a decay exponent 



(oo) 3 - \/5 

V± = — ^ — 



(18) 



in the strong disorder limit. 

Finally, outside the critical point the mapping is re- 
lated to a biased RW, with a finite drift velocity, which 
is proportional to t. From the surviving probability of bi- 
ased RW-s one obtains for the correlation length critical 
exponent [p2|: 



,(oo) 



(19) 



The scaling exponents and relations in Eqs.(p^, (17), 
(|l|) and (|W|) are identical with those of the IRFP of the 
one-dimensional RTIM ||l2|, which is known to control 
the critical behavior of several other random quantum 
spm chains ||l|,|l§ and non-equilibrium phase transitions 
in the presence of quenched disorder p8| . At this point 
our next question is about the region of attraction of 
the IRFP. For the RTIM, where the RW mapping can be 
generalized for weaker disorder, any small amount of ran- 
domness seems to bring the system into the IRFP [ p^ , 
which claim is checked by intensive numerical calcula- 
tions |2^,|2^,^. There are, however, several other mod- 
els (random quantum clock-model, Ashkin- Teller model 



|19| , directed percolation 
romagnetic spin chains |l6[l 



8|, S* = 1 random antifer- 
etc.) where weak disorder 



is not sufficient to bring the system into the IRFP. In 
these cases either the pure systems fixed point stays sta- 
ble against weak disorder perturbations or the competi- 
tion between (quantum) fluctuations and weak quenched 
disorder leads to conventional random scaling behavior. 
For the random percolation problem the latter scenario 
is likely to happen, since the RW mapping can not be 
extended for small disorder. (The transformation law for 
the connected sites, nk/nk-i ~ q ot 1/q, does not hold 
aromid q « 1/2.) We are going to study this issue nu- 
merically by MC simulations in the next Section. 



IV. NUMERICAL RESULTS 

For intermediate strength of disorder we studied the 
percolation by MC simulations. Since the critical proper- 
ties of the problem are related to the connectivity prop- 
erties of clusters for this purpose we implemented the 
standard Hoshen-Koopelman labelling algorithm |^ . To 
decide about the shape of the lattice one should take 
into account the expected anisotropic scaling properties 
of the system, since the scaling functions, as in Eq.(|^) 
depend on the ratio r = L'^ / K, where z is an unknown 
parameter. To overcome this difficulties we used a strip- 
like geometry, when K ^ L, thus r « for all strip 
widths. In practice we had K = 10^, went up to i = 64 



and imposed periodic b.c. in the vertical direction. For 
the distribution of the disorder we used the power-law 
form in Eq. ( p^ , which has already been turned out suc- 
cessful in similar investigations for random quantum spin 
chains | p9[ . Since averaging in the vertical direction in 
Eq.(^ (and also in the horizontal direction for bulk cor- 
relations) is equivalent to a partial average over quenched 
disorder it was enough to consider only a limited number 
(~ 10 — 20) realizations. 




FIG. 4. Estimates for the anisotropy exponent for differ- 
ent strength of disorder. The straight lines connecting the 
points are guide to the eye, for D > Doa ~ 1.2 — 1.5 the 
anisotropy exponent is possibly divergent. In the inset extrap- 
olation of the size-dependent effective anisotropy exponents is 
shown, for D = 0.25, 0.5, 0.75 and 1.0, up to dawn. 

First, we determine the anisotropy exponent, z, by cal- 
culating the probability distribution of clusters in Eq.(||). 
While the decay exponent, r in Eq.(^ has only a weak 
anisotropy dependence, the scaling function R{y) turned 
out to be sensitive of the value of z. As we noticed in the 
numerical calculations R{y) has two different regimes. 
For smaller values of the parameter, y = N/L^ < y* , 
the finite size effects are negligible and the scaling func- 
tion is approximately constant. For y > y* , when the 
largest clusters touch the boundaries, the scaling function 
has a characteristic variation. Measuring the position of 
y* for different widths, L, we obtained a series of effec- 
tive anisotropy exponents, which are then extrapolated 
to L — > oo, as shown in the inset to Fig. |^. This proce- 
dure is repeated for several disorder parameters and the 
extrapolated anisotropy exponents are plotted in Fig. U, 
Unfortunately, with this method we could not go to very 
strong disorder, while the cross-over region can not be 
clearly located for D > 1. However, it is clear from the 
available data that z is monotonically increasing with the 
strength of disorder and it is likely that z will be diver- 
gent for D > Doo ~ 1.2 - 1.5. 

In order to obtain more information about the criti- 
cal behavior of the system we have calculated the bulk 
and the end-to-end average correlation functions at the 
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critical point, as defined in Eqs.(^) and (^), respectively. 
In Fig. I the average bulk correlations, C5^(L), vs. L is 
drawn in a log-log plot. The slope of the curves, which 
is related to the decay exponent, r]±, has a disorder de- 
pendence. 
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FIG. 5. Average bulk correlations vs. the width of 
the strip for different strength of disorder, from D = to 
D = 1.75 in units of 0.25 from up to dawn. The typical error 
is generally smaller than the size of the symbols for small D, 
whereas for larger D it is at most twice of the size of symbols. 
The straight lines are least-square fits. 
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FIG. 6. Bulk (ryx) and surface (r;5_) decay exponents ver- 
sus the strength of disorder. Values at the IRFP, as given in 
Eqs.(p^ and (|l^) are denoted by dashed lines. Two typical 
error bars are also indicated. 

The exponents, calculated in this way together with 
the decay exponent of the end-to-end correlations, 77;^, 
are plotted in Fig. ||. As seen in Fig. |^ both exponents 
are monotonously increasing with the strength of disorder 
and tend to saturate at the respective IRFP values, given 
in Eqs.(|l8|) and (p^. The value of disorder strength, 
where the saturation takes place, within the error of the 
calculation, is the same for the two exponents and it is 



compatible with the estimate, Doo, as calculated from 
the divergence of the dynamical exponent in Fig. ^. 

We can thus conclude that the critical behavior of the 
random percolation process has a weak-to-strong disor- 
der cross-over. For weaker disorder, D < Doo, what we 
call the intermediate disorder regime, the critical behav- 
ior of the system is controlled by a line of conventional 
fixed points. Here the anisotropy exponent is finite, and 
together with the order-parameter exponents, r/± and rj']_, 
monotonously increasing with the strength of disorder. 
In the strong disorder regime, D > Doo, the critical be- 
havior of the system is controlled by the IRFP. Here the 
anisotropy exponent is formally infinity and the other 
critical exponents have no disorder dependence. 



O 

Ln 

CO 



1.2 



1.1 



1.0 



o 



A 



□ 



o 



L=16 
L=24 
L=32 
L=40 



□ 
O 

A 



□ 



□ 



o 



□ 



A 



o 



-0.2 



-0.1 



tL 





1/2 



0.1 



0.2 



FIG. 7. Scaling plot of the bulk correlation function with 
h'± — 2 ai a. disorder strength D = 0.75. 

The non-universal nature of the critical behavior in the 
intermediate disorder regime is possibly connected to the 
presence of a marginal operator, which should have van- 
ishing anomalous dimension, Xe = 0, in the entire disor- 
der range, < D < Dao- In our case the disorder per- 
turbation is connected to the local energy-density opera- 
tor, for which the marginality condition, according to the 
Harris criterion in Eq.([l^) requires the condition (f> = Q, 
thus i/j^ = 2. To verify this scenario we have calculated 
the average bulk correlation function, C\{L,t), outside 
the critical point, at a disorder strength, D — 0.75, which 
is in the middle of the intermediate disorder regime. Ac- 
cording to scaling considerations 



Ci(L,t)=L-'^^C{tL^/''^) 



(20) 



thus from an optimal scaling collapse v± can be de- 
termined. As shown in Fig. |^ the scaling behavior 
of C\{L,t) is compatible with the conjectured value of 
i/j^ = 2 and thus with the marginality condition. 
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V. DISCUSSION 

In this paper bond percolation is studied on the square 
lattice with strictly correlated, layered randomness. The 
phase diagram of the problem as a function of the 
strength of disorder contains two regions. For strong 
enough disorder the critical properties of the model are 
controlled by an IRFP, the properties of which are ex- 
actly known by a RW mapping. For weaker disorder, in 
the intermediate disorder regime the critical behavior is 
found to be controlled by a line of conventional random 
fixed points, where both the anisotropy exponent and the 
order-parameter exponents are disorder dependent. The 
correlation length exponent, however, stays constant at 
its marginal value. 

This type of critical behavior is very similar to that ob- 
tained in a class of random quantum spin chains ||l9|,|l8| . 
This close similarity can be understood by noting the re- 
lation between percolation and the Q — > 1 limit of the Q- 
state ferromagnetic Potts model. With layered random- 
ness the two-dimensional Potts model in the Hamiltonian 
limit p^ ] is equivalent to a quantum Potts chain, with 
random coupHngs, Ji, and transverse fields, hi. ||2^, the 
critical behavior of which can be studied by the SDRG 
method Q. In this procedure the couphngs and trans- 
verse fields are put in descending order and the strongest 
terms are successively decimated out, whereas neighbor- 
ing terms are replaced by renormalized values. Decimat- 
ing the strongest coupling, say J2, yields a new effective 
spin cluster in a renormalized transverse field of strength: 

where hi and /12 are the original transverse fields act- 
ing at the two end-spins of J2. Similarly, if, the spin 
in the strongest transverse field, /i2, is decimated out, 
then a new renormalized coupling is generated between 
remaining spins, which is of the form in Eq. (pi]) , by in- 
terchanging hi <-^- Ji, which is due to duality. 

If the disorder is strong enough, so that the system 
under renormalization is in the attractive region of the 
IRFP, the model specific prefactor 2/Q in Eq.(|2|) does 
not matter and the critical properties are universal. The 
region of strong attraction of the IRFP, however, is lim- 
ited by Q = 2, i.e. for the RTIM. For smaller values of 
Q, like in percolation, when the prefactor in Eq.(^) is 
larger than one, for weak disorder some renormalized cou- 
plings and transverse fields are larger than the decimated 
ones. If this happens frequently, i.e. when the disorder 
is too weak, then the SDRG method is no longer valid 
and the critical behavior of the model is expected to be 
controlled by a conventional random fixed point. This is 
exactly what we obtained by MC simulations. 

We close our paper with two remarks. First, for 
strong enough disorder the critical behavior of both or- 
dinary and directed percolation ||l^ is controlled by the 



same IRFP, thus the original anisotropy between the two 
pure problems does not make any influence about the 
(strongly) random critical behavior. Our second remark 
concerns possible Griffiths effects in the random percola- 
tion problem. Using the analogy with random quantum 
spin chains for strong disorder some dynamical quanti- 
ties of the random percolation problem are singular also 
outside the critical point. For example the susceptibil- 
ity in a uniform field, H diverges as x ~ //-i+i/^ ^ and 
the vertical correlation function decays algebraically as 
C||(/) '-^ , where z' is a finite dynamical exponent, 

which depends on the distance of the critical point. 
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